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Angle resolved photoemission spectroscopy (ARPES) provides a detailed view of the renormal- 
ized band structure and, consequently, is a key to the self-energy and the single-particle Green's 
function. Here we summarize the ARPES data accumulated over the whole Brillouin zone for the 
optimally doped Bi2Sr2CaCu20s-5 into a parametric model of the Green's function, which we use 
for calculating the itinerant component of the dynamic spin susceptibility in absolute units with 
many-body effects taken into account. By comparison with inelastic neutron scattering (INS) data 
we show that the itinerant component of the spin response can account for the integral intensity of 
the experimental INS spectrum. Taking into account the bi-layer splitting, we explain the magnetic 
resonances in the acoustic (odd) and optic (even) INS channels. 

PACS numbers: 74.72.-h 74.72.Hs 74.25.Ha 74.20.-z 79.60.-i 



The origin of the magnetic resonance structure observed 
in the superconducting (SC) state of YBa 2 Cu30 6+( 5 
(YBCO) 0, i i |, [|, Bi 2 Sr 2 CaCu 2 0= 5 (BSCCO) 
@i 0, S] , and other families of cuprates Q is one of the 
most controversial topics in today's high-T c superconduc- 
tor (HTSC) physics. Existing theories waver between 
the itinerant magnetism resulting from the fermiology 
[lol . [ill . [H, [HI, EH and the local spins pictures (such as 
static and fluctuating "stripes" , coupled spin ladders, or 



and compare the resulting spectrum calculated for op- 
timally doped BSCCO with the available INS measure- 
ments on both BSCCO and YBCO. 

We start with establishing the relation between the 
quasiparticle Green's function and INS response. The 
normal-state Lindhard function is related to the quasi- 
particle Green's function via the following summation 
over Matsubara frequencies [TU, [l9[ : 



spiral spin phase models) 

both approaches can qualitatively reproduce the main 
features of the magnetic spectra in the neighborhood 
of the optimal doping. It is a long standing question, 
which one of these two components (itinerant or local- 
ized) predominantly forms the integral intensity and the 
momentum-dependence of the magnetic resonances. It is 
therefore essential to estimate their contribution quanti- 
tatively, carefully taking into account all the information 
about the electronic structure available from experiment. 
However, such a comparison, which could shed light on 
the dilemma, is complicated, as it requires high-quality 
INS data and the extensive knowledge of the electronic 
structure for the same family of cuprates. On the other 
hand, APRES data for YBCO compounds, for which the 
best INS spectra arc available, are complicated by the 
surface effects 17], while for BSCCO, most easily mea- 
sured by surface-sensitive techniques such as ARPES, the 
INS measurements show much lower resolution due to 
small crystal sizes. 

Here we propose a way to estimate the dynamic spin 
susceptibility in the odd (o) and even (e) channels within 
the random phase approximation (RPA) from the single- 
particle spectral function, including many-body effects, 
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Besides the bare Green's function, equation Jl} also holds 
for the renormalizcd one. It can be rewritten as a double 
integral along the real energy axis [2(| [2l[ : 
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where C tj (k, e,v) = Im G i ( k > e ) Im G j {k + Q,v)dk 
is the cross-correlation of the constant-energy cuts of the 
spectral function over the Brillouin zone (BZ), H/(e) = 
l/(e e / fcsT + 1) is the Fermi function, and indices i and j 
numerate the bonding and antibonding bands. The fac- 
tors Cij (k, e, v) can be efficiently calculated in the Fourier 
domain by means of the cross-correlation theorem pjj . 

In the SC state, the anomalous Green's function F 
additionally contributes to xo [3 : 

dj (k, e,u) = [Im Gi (k, e) Im G 3 (k + Q, u) 

+ ImF i (k,e)ImF j (k + Q,i/)] dk (3) 

Although ImF is not directly measured by ARPES, it 
can be still accounted for, as we will subsequently show. 
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Fig. 1 (color online). Comparison of the model with experimental ARPES spectra of OD Bi-2212 at 30 K taken with 50 eV (a, b) 
and 38 eV (c) photon energies. The model spectra are smoothed with a gaussian to account for 20 meV energy resolution and 
0.025 A -1 angular resolution, a) Spectra at the (7r, 0) point with the corresponding energy distribution curves (below) taken 
along the dashed lines, b) Nodal spectra along the (n, n) direction, c) Comparison of the experimental and model spectra 
taken at an intermediate position in fc-space to check the validity of the interpolation of the self-energy between the nodal and 
antinodal directions, d) The Kramers-Kronig consistent real and imaginary parts of the nodal and antinodal self- energies, e) 
Positions of the cuts a)-c) in fc-space. 



After one knows the Lindhard function xo (frequently 
referred to as the bare spin susceptibility) , one can fi nally 
obtain from RPA the dynamic spin susceptibility x 
the imaginary part of which is directly proportional to 
the measured INS intensity 

X °-°(Q, 0) = Xo'°(Q, O)/ [1 - J° c xS' e (Q, O)] (4) 

The coefficient Jq ° in the denominator of (0| describes 
the effective Hubbard interaction. In our calculations we 
employed the model for Jg' e 
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-J||(cosQ x + cos Q v ) ± Jj_, (5) 



where the first term accounts for the Q-dependence due 
to the in-plane nearest-neighbor superexchange, and the 
second term arises from the out-of-plane exchange inter- 
action. 

Thus, knowing the single-particle Green's function 
leads us to a comparison of ARPES results with the 
INS data. The previous calculations based on this idea 
were performed for the bare band struc- 
ture only, disregarding the rcnormalization effects, which 
makes the conclusions based on comparison with the INS 
data rather uncertain. The recent work by U. Chatterjee 
et al. [24} is the only available paper that includes the 
many-body effects from experimental data (in a proce- 
dure different from ours), but it does not account for the 
bi-layer splitting (necessary for reproducing the odd and 
even INS channels), provides the results in arbitrary units 
only, rather than on an absolute scale, and gives only an 
estimate for the anomalous contribution to x- So we will 
address these issues in more detail below. 



At first, we will introduce an analytical model that can 
reproduce the ARPES measurements within a wide en- 
ergy range and all over the BZ. As in a single experiment 
it is practically impossible to obtain a complete data set 
of ARPES spectra, such a model allows making use of all 
the available data measured from a particular sample and 
calculating the full 3D data set afterwards. In such a way 
the effect of matrix elements and experimental resolution 
is also excluded. 

The measured ARPES intensity is basically propor- 
tional to the imaginary part of the Green's function (al- 
though it is affected by matrix elements, experimental 
resolution, and other factors [IBf). The latter can be 
obtained if one knows the self-energy, extracted from the 
ARPES data in a routine self-consistent Kramers-Kronig 
procedure (26[. 

We employed a model of the Green's function based on 
the bare electron dispersion studied in [2(| and a model 
for the imaginary part of the self-energy E" = E"j + E[J os , 
where E"j = aoj 2 is the Fermi- liquid component of the 
scattering rate that originates from the electron-electron 
interactions, and E(^ os models the coupling to a bosonic 
mode [13] ■ In the (tt, 7r) (nodal) direction we modeled 
E£ os by a step-like function E£ Q = \ & [l+tanh(=^)] 
of width (5w n , height f3 n and energy £l n , while in the (w, 0) 
(antinodal) direction we accounted for the peak in the 
self-energy due to the pile-up in the density of states 
at the gap energy: E," = — /3 a Rc = , 

6 f S>J bos ^(w-i5w a ) 2 -(Ao + I2 a ) 2 ' 

where Ao is the SC gap, f2 a is the mode energy, and 
5u> a is the broadening parameter (see [28} and references 
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Fig. 2 (color online), a) Energy dependence of the real and 
imaginary parts of the Lindhard function xo at the (n, n) 
point for the odd and even channels, b) Contributions of 
the normal (thin solid curves) and anomalous (dashed curves) 
components to the real and imaginary parts of Xo dd in the SC 
state. The sum of two components is shown as thicker curves. 
In our calculations we used the F value in ([2]) of 5 meV, which 
could introduce insignificant additional broadening of xo as 
compared to the bare band calculations. The energy integra- 
tion range in ([2]) was chosen to be ±0.25 eV. 



therein) . The real part of the self-energy £' was then de- 
rived by the Kramers-Kronig procedure and the Green's 
function was calculated according to [291]: 
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where A(k) is the SC d-wave gap changing from zero 
along the BZ diagonals to the maximal value of ±Ao 
along the antinodal directions. Self-energy parameters 
were specified independently for the nodal and anti-nodal 
parts of the spectra, with a d-wave interpolation between 
these two directions: S"(k,w) = E^(w) + j^'M - 
S^w)] (cosk x — cosfcy) 2 . We also assumed the particle- 
hole symmetry in To achieve the best reproduction 
of the experimental data, all the free parameters were 
adjusted during comparison with a set of ARPES spectra 
of Bi-2212 to achieve the best correspondence (Fig.[T]). 
The best-fit parameters of the model are listed in the 
following table: 

a = 3.0 eV -1 f3 n = 30meV /3 a = 200 meV 5u n = lOmeV 
<5u; a = 0.08 A Q„ = 60 meV H a = 42 meV A = 35 meV 

We would like to stress here that such a simple self- 
energy model that includes coupling only to a single 
bosonic mode can accurately reproduce the state of the 
art ARPES spectra of BSCCO, as we have just shown. 

The described model has multiple advantages for nu- 
merical calculations over raw ARPES data. Only in such 
a way one can completely separate the bonding and an- 
tibonding bands, which is impossible to achieve in the 
experiment, and reveal the nature of the odd and even 
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energy S(k, uj) and the pairing vertex A(k) 

present in the model, it becomes possible to calculate the 
anomalous Green's function F(k,uj) [30l ] : 
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[w-S(k,w)] - A 2 (k) 

Besides the already mentioned absence of matrix element 
effects and experimental resolution, the formulae (O and 
([7]) also allow to obtain both real and imaginary parts 
of the Green's functions for all k and lu values including 
those above the Fermi level. It automatically implies the 
particle-hole symmetry (ekj-k = — ekf+k) in the vicinity 
of the Fermi level, which in case of the raw data would 
require a complicated symmctrization procedure based 
on Fermi surface fitting, being a source of additional er- 
rors. Finally, it provides the Green's function in absolute 
units, allowing for quantitative comparison with other ex- 
periments and theory, even though the spectral function 
originally measured by ARPES lacks the absolute inten- 
sity scale. Thereupon, we find the proposed analytical 
expressions to be better estimates for the self-energy and 
both Green's functions and therefore helpful in calcula- 
tions where comparison to the experimentally measured 
spectral function is desirable. 

Now we will apply the described model to calculating 
the dynamic spin susceptibility. Starting from the model 
data set built for optimally doped BSCCO at 30 K, with 
the maximal SC gap of 35 meV, we have calculated the 
Lindhard function (Eq.[5]) in the energy range of ±0.25 eV 
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Fig. 3 (color online), a) k- integrated \" = Im(x) in the 
odd (solid curve) and even (dashed curve) channels, b) k- 
dependence of the resonance energies in odd (*) and even (a) 
channels along the high-symmetry directions (0, 0) - (77, 77) - 
(77,0). The dashed lines mark the onset of the particle-hole 
continuum (position of the "step" in Xo)- Second row: Con- 
stant energy cuts of x" in the odd channel below the resonance 
(c), at the resonance energy (d), and above the resonance (e). 
The center of each BZ image corresponds to the (77, 77) point. 
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in the whole BZ for the odd and even channels of the spin 
response (see Fig. [2^,). To demonstrate that the contri- 
bution of the anomalous Green's function is not negligi- 
bly small, in Fig. [2b we show separately the normal and 
anomalous components of Xo dd - 

After that we calculated x (Eq.EJ by adjusting the J\\ 
and Jj_ parameters to obtain correct resonance energies 
at (it, 7t) in the odd and even channels (42 and 55 meV 
respectively), as seen by INS in BSCCO @, 0, 0- Thc 
resulting x°' ° are qualitatively similar to those obtained 
for the bare Green's function 14[ . The intensity of the 



resonance in the even channel is approximately two times 
lower than in the odd channel, which agrees with the ex- 
perimental data 0, [||]. On the other hand, for Jj_ = 
the splitting between odd and even resonances does not 
exceed 5-6 meV, which is two times less than the experi- 
mental value. This means that the out-of-plane exchange 
interaction (in our case Jl/Jii ~ 0.09) is significant for 
the splitting and the difference in xo alone between the 
two channels cannot fully account for the effect. 

In Fig.[3K we show both resonances, momentum- 
integrated all over the BZ. Here we would like to draw 
the reader's attention to the absolute intensities of the 
resonances. A good estimate for the integral intensity 
in this case is the product of the peak amplitude and 
the full width at half maximum, which for the odd reso- 
nance results in 0.12 /i^/f.u. in our case. This is in good 
agreement with the corresponding intensity in latest ex- 
perimental spectra on YBCO (~0.11 fi%/f.u.) [||. 

As for the momentum dependence of x> Fig.[3j3 shows 
the dispersions of incommensurate resonance peaks in 
both channels along the high-symmetry directions, cal- 
culated from the Green's function model with the self- 
energy derived from the ARPES data. We see the W- 
shaped dispersion similar to that seen by INS on YBCO 
0, EH and to the one calculat ed p reviously by RPA for 



the bare Green's function [13|, [14(. At (tt,tt) both reso- 
nances are well below the onset of the particle-hole con- 
tinuum at ~ 65 meV (dashed line), which also agrees with 
previous observations 0, EH, EH • At higher energies mag- 
netic excitations are ovcrdamped, so the upper branch of 
the "hourglass" near the resonance at (n, it) suggested by 
some INS measurements 0, 0, EH is too week to be ob- 
served in the itinerant part of x and is either not present 
in BSCCO or should originate from the localized spins. 

In Fig. [3] we additionally show three constant-energy 
cuts of x m the odd channel below the resonance, at 
the resonance energy, and above the resonance. As one 
can see, besides the main resonance at (it, it) the calcu- 
lated x reproduces an additional incommensurate reso- 
nance structure, qualitatively similar to that observed in 
INS experiments [2| . Below the resonance the intensity is 
concentrated along the (k,0) and (0, k) directions, while 
above the resonance it prevails along the diagonal direc- 
tions (k, ±fc). 

In this work we have demonstrated the basic relation- 



ship between the ARPES and INS data. The comparison 
supports the idea that the magnetic response below T c 
(or at least its major constituent) can be explained by 
the itinerant magnetism. Namely, the itinerant compo- 
nent of Xi at least near optimal doping, has enough in- 
tensity to account for the experimentally observed mag- 
netic resonance both in the acoustic and optic INS chan- 
nels. The energy difference between the acoustic and op- 
tic resonances seen in the experiments on both BSCCO 
and YBCO, cannot be explained purely by the differ- 
ence in xo between the two channels, but requires the 
out-of-plane exchange interaction to be additionally con- 
sidered. In this latter case the experimental intensity ra- 
tio of the two resonances agrees very well with our RPA 
results. Also the calculated incommensurate resonance 
structure is similar to that observed in the INS exper- 
iment. Such quantitative comparison becomes possible 
only if the many-body effects and bi-layer splitting are 
accurately accounted for. A possible way to do that is to 
use the analytical expressions for the normal and anoma- 
lous Green's functions proposed in this paper. We point 
out that such method is universal and can be applied 
also to other systems with electronic structure describ- 
able within the self-energy approach. 
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